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Abstract 

A systematic calculation of nuclear matter is performed which includes the long-range 
correlations between nucleons arising from one- and two-pion exchange. Three-body effects 
from 27r-exchange with excitations of virtual A(1232)-isobars are also taken into account 
^ : in our diagrammatic calculation of the energy per particle E{kf). In order to eliminate 

possible high-momentum components from the interactions we introduce at each pion- 
■ baryon vertex a form factor of monopole type. The empirical nuclear matter saturation 

point, pq ~ 0.16fm^^, ~ — 16MeV, is well reproduced with a monopole mass of 
A ~ 47r/,r ~ 1.16 GeV. As in the recent approach based on the universal low-momentum 
A^iV-potential T^ow-k; the inclusion of three-body effects is crucial in order to achieve 
saturation of nuclear matter. We demonstrate that the dependence of the pion-exchange 
contributions to E{kf) on the "resolution" scale A can be compensated over a wide range 
of A by counterterms with two "running" contact-couplings. As a further application we 
study the in-medium chiral condensate {qq){p) beyond the linear density approximation. 
\^ I For p < l.bpo we find small corrections from the derivative dE{kf)/dmT^, which are stable 

I against variations of the monopole regulator mass A. 

O ■ PACS: 12.38.Bx, 21.65.+f 

Keywords: Nuclear matter equation of state; one- and two-pion exchange with pion-nucleon 
monopole form factors; in-medium chiral condensate 

1 Introduction and summary 

;_j ■ In recent years a novel approach to the nuclear matter problem based on effective field theory 
has emerged. Its key element is a separation of long- and short-distance dynamics and an order- 
ing scheme in powers of small momenta. At nuclear matter saturation density po — 0.16 fm~^ 
the Fermi momentum kfQ and the pion mass are comparable scales {kfo c:^ 2m^). Therefore 
pions must be included as explicit degrees of freedom in the description of the nuclear many- 
body dynamics. The contributions to the energy per particle, E{kf), of isospin-symmetric 
(spin-saturated) nuclear matter as they originate from chiral pion-nucleon dynamics have been 
computed up to three- loop order in refs. [1,2]. Both calculations are able to reproduce the empir- 
ical saturation point by adjusting one single parameter (either a contact-coupling go + gi — 3.23 
[1] or a momentum cut-off A ^ 0.65 GeV [2]) related to unresolved short- distance dynamics. 
The mechanism for saturation in these approaches is mainly a repulsive contribution to the 
energy per particle generated by Pauli-blocking in second order (iterated) one-pion exchange. 
As outlined in ref . [2] this mechanism becomes particularly transparent by taking the chiral limit 
= 0. In that case the interaction contributions to the energy per particle are completely 
summarized by an attractive fc^-term and a repulsive kj-term where the parameter-free predic- 
tion for the coefficient of the latter is very close to the one extracted from a realistic nuclear 
matter equation of state. 

^Work supported in part by BMBF and GSI. 
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This chiral approach to nuclear matter has been extended and improved in rcf. [3] by includ- 
ing systematically the effects from 27r-cxchange with excitation of virtual A(1232)-isobars. The 
physical motivation for such an extension is threefold. First, the spin-isospin-3/2 A(1232)- 
resonance is the most prominent feature of low-energy vrA/'-seattering. Secondly, it is well 
known that the 27r-exchange between nucleons with excitation of virtual A-isobars generates 
the medium- and long-range components of the isoscalar central iViV-attraction [4]. In phe- 
nomenological one-boson exchange models this vrA^A-dynamics is often simulated by a scalar 
"cr" -meson exchange. Thirdly, the delta- nucleon mass-splitting A = 293 Me V is of the same 
size as the Fermi momentum /c/o — Sm^r at nuclear matter saturation density. Therefore pions 
and A-isobars should both be treated as exphcit degrees of freedom. 

It has been found in rcf. [3] that the inclusion of the chiral vriVA-dynamics significantly 
improves e.g. the momentum-dependence of the (real) single-particle potential U{p,kf ) and 
the isospin properties (as revealed in the density-dependent asymmetry energy A{kf) and the 
neutron matter equation of state En{kn))- However, there remain some open questions in 
these perturbative calculations of nuclear matter. The effective short-range terms are adjusted 
directly to nuclear matter bulk properties and thus seem unrelated to those relevant for free 
A^A^-scattering. Also, no deeper justification for a perturbative treatment of nuclear matter 
(besides of being very successful) could be given. 

Important progress in this direction has come recently from the work of Bogner et al. [5] 
based on the universal low-momentum A'"A'"-potential VJow-k- This potential operates by con- 
struction only between the low-momentum nucleon states, |p| < 0.4 GeV, where it is truly 
determined by elastic A^iV-scattering data. It has been demonstrated in ref. [5] that, due to the 
absence of a model-dependent short-range repulsive core in the potential Viow-k, its second and 
higher order iterations in nuclear matter (i.e. the successive terms in the Brueckner ladder- 
series) turn out to be small. The primary reason for this (unconventional) feature is the Pauli 
exclusion principle: with increasing Fermi momentum kf Pauli-blocking reduces progressively 
the available low-momentum phase space wherein Viow~k acts. Under such conditions a per- 
turbative treatment of the two-body interaction in nuclear matter turns out to be justified. It 
has also been demonstrated in ref. [5] that a saturation of nuclear matter cannot be achieved 
from the (low-momentum) two-body interaction alone (see also ref. [6] for earher Hartree-Fock 
calculations with VJow-k exhibiting the same feature). Repulsive three-nucleon interactions (in 
particular the long-range ones induced by 27r-exchange) are crucial in order to stabilize the 
nuclear many-body system against collapse. As emphasized in ref. [5] in an effective low-energy 
theory (with a "spatial resolution" of A = (0.4 GeV) ~^ ~ 0.5 fm as it is inherent to the Viow-k 
potential) many-body forces are inevitable. Taking into account furthermore the dominant 
second order contributions from VJow-k and the three-nucleon force (adjusted partly to triton 
and "^He binding energies) a saturation minimum of Eq ^ — 12MeV at the correct equilibrium 
density kfQ ~ 1.35 fm^^ has been obtained [5]. Presently, there is ongoing activity to further 
improve these calculations in order to close the remaining small gap to the empirical binding 
energy per particle, Eq ^ — 16MeV. 

The purpose of the present paper is to reanalyze the chiral approach to nuclear matter 
from this new perspective and to establish connections to the l^ow-k approach (which actually 
starts from phase-shift equivalent A^A^-potentials). As already mentioned, the second-order 
tensor force from iterated one-pion exchange plays an essential role in order to obtain binding 
and saturation of nuclear matter in the (perturbative) chiral approaches studied so far [1, 2, 
3]. On the other hand its strong short-distance components are one of the sources which 
render the Brueckner ladder-series in nuclear matter non-perturbative (i.e. not convergent) [5]. 
It is therefore advisable to eliminate the high-momentum components from the pion-induced 
A^iV-interactions. For that purpose we follow the phenomenology of the one-boson exchange 
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Figure 1: Two- loop one-pion exchange Fock diagram and three- loop iterated one-pion exchange 
Hartree and Fock diagrams. The combinatoric factors are 1/2, 1/4 and 1/4, in the order shown. 
Their isospin factors for symmetric nuclear matter are 6, 12 and —6, respectively. 



potentials and introduce at each pion-nucleon interaction vertex a "form factor" of monopole 
type: 



A2 
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with A the monopole mass. Here, < denotes the squared four-momentum transfer carried 
by the virtual pion. Typical values of the monopole mass A employed in one-boson exchange 
potentials lie in the range: 1.0 GeV < A < 1.7 GeV [7, 8]. With the form factor F{q^) included 
as a regulator on high-momentum components the Ivr- and 27r-exchange contributions to the 
energy per particle E{kf) depend on the monopole mass A. As a first orientation we find that 
the "natural choice" A ~ 4:71 f^^ ~ 1.16 GeV (the chiral symmetry breaking scale) reproduces 
correctly the empirical nuclear matter saturation point: po = 2fcjQ/37r^ ~ 0.16 fm~^, E{kfQ) ~ 
— 16MeV. As in the V^ow-k approach of ref.[5] the inclusion of the pion-induced three-body 
terms (with and without virtual A-isobar excitation) is essential in order to achieve saturation 
of nuclear matter. 

However, a physical low-energy quantity such as the nuclear matter equation of state E{kf) 
should not depend on a parameter A which merely scans the "spatial resolution" of the pion- 
baryon interactions involved. We demonstrate that this (unphysical) A-dependence of the 
pion-exchange contributions is, over a large range of monopole masses and densities, almost 
perfectly counterbalanced by the fc^- and fcj-terms related to two "running" short- distance 
contact-couplings ,835 (A). 

Having reached a description of nuclear matter in terms of explicit pion-exchange dynamics, 
which is stable against variations of the short-distance scale A, we can study the in-medium chiral 
condensate {qq) (p) beyond the linear density approximation. These corrections are obtained by 
differentiating the calculated energy density of nuclear matter pE{kf) with respect to the pion 
mass (or equivalently, the light quark mass m^). Below the saturation density p < 0.16 fm~^ we 
find very small corrections. At higher densities a tendency counteracting chiral restoration sets 
in. Moreover, there is little dependence of the derivative dE{kf)/dmT^ on the monopole mass A 
(assuming the short- distance contact-couplings ^3^5 (A) to be independent of the quark mass). 

Our paper is organized as follows: In section 2 we present first the analytical expressions 
for the one- and two-pion exchange contributions to the energy per particle E{kf) including the 
monopole form factors. Then we discuss in section 3 the results for the nuclear matter equation 
of state together with the compensating short- distance contact-couphngs 535(A). Section 4 is 
devoted to the in-medium chiral quark condensate {qq){p)- Section 5 ends finally with some 
concluding remarks and an outlook. 
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2 Diagrammatic calculation of the energy per particle 



The first contribution to the energy per particle E{kf ) is the kinetic energy of relativistic Fermi 
gas of nucleons (expanded in powers of 1/Mjv): 

^/K N(kin) _ , ^/ (2) 



Here, kf denotes the Fermi momentum, related to the nucleon density in the usual way, p 
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2ky3TT^, and = 939 MeV stands for the (average) nucleon mass. At the densities of interest. 



p < 0.5 fm~^, the last term in eq.(2) is already negligibly small (< 0.1 MeV). Next, we enumerate 
the interaction contributions from one- and two-pion exchange. For each diagram we present 
only the final result omitting all technical details related to solving the loop and Fermi sphere 
integrals in the presence of the monopole form factors. 

2.1 One-pion exchange Fock term 

The contribution of the left l7r-exchange Fock diagram in Fig. 1 (including the relativistic 
correction) reads: 

- |g{^:^^^ + arctan(2.) + i^arctan(2rK) 

1 + 12«2 3r2 - 2 + 12rV(2r2 - 1) ^ ^ 



ln(l + 4u^) H ^ ln(l + 4r^u^ 



32^3 ^ ^ 32rS3 



ml [r^ - + 3 3m 



40Mjf 



2r^u r 



^ (3r^ - 7r^ + 4) - (5 + 12^^) arctan(2M) 



+-^(25r^ - 15 - 12rV + 36rV) arctan(2rw) - ^ ln(l + Au'^) 

4r^ — 3 
+— ^^ln(l + 4rV) 



(3) 



For notational simplicity we have introduced the dimensionless variable u = kf/rriT^ and the ratio 
r — m^/A. The other occurring parameters are: qa = 1-3 (the nucleon axial- vector coupling 
constant), = 135 MeV (the neutral pion mass), and /,r = 92.4 MeV (the weak pion-decay 
constant). In the limit r — > one recovers the result eq.(6) in ref.[2] without the monopole form 
factor. 

2.2 Iterated one-pion exchange 

The Hartree diagram of iterated Ivr-exchange (middle diagram in Fig. 1) with two medium 
insertions^ leads to the following analytical expression: 

F,/, Nm2i 3ylMjvm^ f 20(3 + r^) + 32(1 + r2)M2 

^^^^ " 5(87r)V I ^[arctan(2«)-arctan(2r«) 

1/1 5 11 41 4 \ / 15 41 75r 
+ - - + -r^-:r^ + — + - + u 29 - 



u 



2 16r^ 8r3 16r 1 + r/ V Ar 8 1 + r 



^Medium insertion is technical notation for the difference between the nucleon propagator in the medium and 
in the vacuum (see sect. 2 in ref.[2]). To each medium insertion belongs an integration over a Fermi sphere of 
radius kf. 
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r2-9-20(7 + r2)M2 



+ 



8(1 -r2)M3 
27r2 - 5 - 63r^ + lOSr^ 



ln(l + W) 



(4) 



64r7(l _ r2)M3 



ln(l + 4r^'u^) 



The corresponding exchange term comes from the Fock diagram of iterated Ivr-exchange (right 
diagram in Fig. 1) with two medium insertions. While the loop integral including the two 
squared monopole form factors can be solved there remains a single non-elementary integral 
from the integration over the Fermi spheres. We get: 



E{kf) 



(F2) _ 3fflMjvm^ 



ru 

/ dx xiu — xYi2u + x) 
Jo 



1 + Sx^ + 8x^ 
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X 
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arctan ■ 
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2rx 



1 + r 



2 l + r^ + 4rV(l + r2 + 2rV)"' 



1 — r 
arctan(rx) 



4r2(l + 2r2a;2)3 
1 + + 4r*x^{l + r^ + 2r^x^) 1 + + ^r^x^ 



+x(l 



4r2(l + 2r2x2)3 
1 + 



+ 



(1 + r2 + 4r2x2)2, 
(1 - rf + Ar^x^ 



arctan(2ra;) 



1 + r2 + 4r2x2 (1 + r2 + 4r2x2)2 
4^^2(^2 _ 2) + 5,^ _ 9,^-1 + rx^{\ + ,^2) 



4r2 



16(l + r2x2)2 (l + 2r2x2)2 Jj- 

In our way of organizing the many-body calculation, the Pauli-blo eking effects arc represented 
by diagrams with three medium insertions [2, 3]. The contribution of the Hartree diagram with 
three medium insertions (including the fourth power of a monopole form factor) reads: 



(H3) 



dxx" 



dy 



2„.2\ 



(471^)4^3 J, 

|t^[ln(l + rV)-ln(l + s2) 



X 



2uxy + {u — X y 

2s^ 



In 



u + xy 
u — xy 



+ 



■ [6 - + Ibr^s^ 



1 + S2 

8rh^ + 7rS^ + r^s^ 



3rS6 



3(l + r2s2)3r - ■ — " ^. ....... . ^. ^jj, (6) 

with the abbreviation s — xy + \Ju^ — x^ ^ x'^y'^. On the other hand one gets from the right 
Fock diagram in Fig. 1 with three medium insertions (including two squared monopole form 
factors) : 



E{kf) 



(F3) 
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8 



'^^yze{y' + z^-l) 



2„2\ 



-ln(l + s^) +ln(l + r^s 



1 

2\j-2 



lyz^y'-^ + z'^ 



1 + r'^s'^ 



-ln(l + t^) + hi(l + r^t^) 



1 + rH'^ 

with another abbreviation t — xz -\- ^/u^ — x"^ -\- x^z"^ and the auxiliary function: 



(7) 



u 



1 



G = (l-r^)- [l + (w + a;)^][l 



4x 



xf\ In 



"^4x 



2r-^ - + {u' - xy + 2u' + 2x 



In 



\ + {u + xf 

1 + (m - ,X)2 

1 + r^(M + x)^ 
1 + r2(M — x)2 



(8) 
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Note that the contributions in eqs.(4-7) carry the large scale enhancement factor M^v- It stems 
from the energy denominator of these iterated diagrams which is proportional to the difference 
of small nucleon kinetic energies. 



2.3 Irreducible two-pion exchange: Spectral representation 

Next, we come to the irreducible 27r-exchange contributions (with no, single, and double A(1232)- 
isobar excitation). Their direct evaluation requires the one-loop A^A^-scattering T- matrix with a 
monopole form factor attached to each pion-baryon vertex. The essential long-distance informa- 
tion about these 27r-exchange one-loop diagrams is already contained in the spectral functions 
(or imaginary parts). The effect of the monopole form factors^ on the spectral functions can be 
easily inferred from the following partial fraction decomposition: 



A--m!\2 1 1 A2-//, 
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— \ A? — J — q^ A? — q^ [A? — q^)' 



(9) 



It instructs us that a pion-exchange with monopole form factors attached to the vertices is 
identical to a point-like pion-exchange minus additional exchanges involving a " heavy particle" 
of mass A. By unitarity this implies that the spectral function of a one-loop 27r-exchange 
diagram with monopole form factors has thresholds at = 2m7r, /i = A + and ji = 2A. 
At the higher thresholds the effective chiral Lagrangian is, however, no more applicable. In 
this case a meaningful and physically reasonable way to include effects from the monopole form 
factor is to cut off the spectral integral a,t /j, — A + m-;^. Putting all pieces together [3] this 
procedure leads to the following 27r-exchange two-body Fock term: 

where ImV^, ImVFc) IniVr and ImliV are the spectral functions of the isoscalar and isovector 
central and tensor A^A^-amplitudes, respectively. Explicit expressions of these imaginary parts 
for the contributions of the triangle diagram with single A-excitation and the box diagrams with 
single and double A-excitation can be easily constructed from the analytical formulas given in 
section 3 of ref. [4]. Note that the /i- and /cj-dependcnt weighting function in eq.(lO) involves 
two subtractions. This way it is guaranteed that the spectral integral at low and moderate 
densities is dominated by low invariant TTTT-masses, 2m7r < < 1 GeV. 

The contributions to the energy per particle from irreducible 27r-exchange (with only nucleon 
intermediate states) can also be cast into the form eq.(lO). The corresponding non- vanishing 
spectral functions read: 



\ jj? - Ami 

ImWr = — 

37r/.(4A)4 



4m^(l + Agl - hg\) + ^,\2Zg\ - Wa - 1) + ^^^^"^^ 



jj? — Ami 



(11) 



^We choose the same monopole form factor for each pion-baryon vertex. 
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Figure 2: Hartree and Fock three-body diagrams related to 27r-exchange with single virtual 
A-isobar excitation. For symmetric nuclear matter the isospin factors are 8, and 8, in the 
order shown. The combinatoric factor is 1 for each diagram. 

2.4 Three-body terms with intermediate A(1232)-isobars 

Finally, we come to the additional three-body terms which arise from Pauh-blocking of interme- 
diate nucleon states. The corresponding closed Hartree and Fock 27r-exchange diagrams with 
single virtual A-isobar excitation are shown in Fig. 2. Their isospin factors are 8, and 8, in 
the order shown. For the (left) three-loop Hartree diagram the integral over the product of 
three Fermi spheres of radius kf can be solved in presence of the monopole form factors and its 
contribution to the energy per particle reads: 




4rV + 'U^(2r2 - 1 + 15r^) (5 + 3r^)u^ , , 

+ 4r^(l + W) + "^^'""^'^^ 

' ^ ^arctan(2rM) }, (13) 



8r3(l - r2) 

with u — kf/rriT^ and r = m7r/A. The nonrelativistic A-propagator shows up in this expres- 
sion merely via the (reciprocal) mass-splitting A = 293 MeV. Furthermore, we have inserted 
in eq.(13) already the empirically well-satisfied relation Qj^na = 'ig-wN/V^ together with the 
Goldberger-Treiman relation (?^jv = Qa^n/ f-w = 13.2. Finally, the contribution of the right 
three-body Fock diagram in Fig. 2, with two squared monopole form factors included, is given 
by: 

^(^^)'""'-4Aigvr<'K2Gi+Gi), (14, 

where we have introduced the two auxiliary functions: 

Gs = 4:X [ arctan(M + x) + arctan(M — x)] + (x^ — — 1) In — 

1 + [u — xy 

2x 

-\ — ^(1 — 3r^){ arctan[r(xi -|- x)] + arctan[r(ii — x)]} 

+{u' -x' + 2r-' - r-') In , (15) 

1 -I- M — xV 



-^-^ [1 + (m + xf] [1 + (m - xf] In ^ ^ 



1 + (w - xf 
x^ + {l- 2,u^)x^ + 2m^ + 2,u^ + r"^ 
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Figure 3: Energy per particle E{kf) of isospin-symmetric nuclear matter versus the nucleon 
density p — 2A;|/37r^. In the dashed line all three-body terms have been omitted. 



u 



-2r-' - —{3 + u 



(1 -2r2) + 



2 - 3r2 
r^x"^ 



In 



1 + r^{u + xY 
1 + r'^{u — x)' 



(16) 



Evidently, the three-body Fock term in eq.(14) is attractive and reduces, by about 1/3, the 
stronger repulsive Hartree term in eq.(13). 



3 Results for the nuclear matter equation of state 

We are now in the position to present numerical results for the nuclear matter equation of 
state. The interaction contributions to the energy per particle E{kf) written down in section 2 
depend on a free parameter, the monopole mass A. A first possible option is to adjust it to (one 
coordinate of) the empirical saturation point. Imposing a minimum of the saturation curve 
E{kf) at kfo — 262 McV (corresponding to an equilibrium density of po = 0.158 fm~^) fixes 
the monopole mass to Aq = 1.145 GcV. This value is surprisingly close to the chiral symmetry 
breaking scale A^ = iTrf-,,- ~ 1.16 GeV. It is also compatible with the typical monopole masses 
employed in one-boson exchange potentials [7]. With that fixed Ao-value the energy per particle 
at the saturation minimum comes out as Eq — —16.7 MeV, in good agreement with the empirical 
value Eq = (— 16 ± 1) MeV. This quite a non-trivial result. 

The decomposition of Eq into contribution from the kinetic energy, Ivr- and 27r-exchange, 
as well as two- and three-body terms is also interesting: = (21.64 -|- 13.63 — 72.86 — 6.41 + 
19.06 — 1.46-1- 4.63 -I- 7.72 — 2.67) MeV, where the nine entries correspond to the terms in eqs.(2- 
7,10,13,14), in that order. The full line in Fig. 3 shows the resulting nuclear matter equation of 
state as a function of the miclcon density p = 2A:']/37r^ up to about 3po — 0.5 fm~^. The curvature 
at its minimum translates into a nuclear matter compressibility of K = kj^f^E'^kfo) = 292 McV. 
It is about 10% larger than the recent extrapolation from giant monopole resonances of heavy 
nuclei, which gave the value K = (260 ± 10) MeV [9]. 
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Figure 4: The "running" contact-couplings -B3(A) and -B5(A) as a function of the resolution 
scale A. 



The inclusion of the 27r-exchange three-body terms eqs.(6,7,13,14) is crucial in order to 
achieve realistic saturation of nuclear matter in our calculation. This feature is in accordance 
with the findings in the work of ref.[5] based on the universal low-momentum A^A^-potential 
Mow-k- It is also compatible with the results of sophisticated many-body calculations by the 
Urbana group [10] . The dashed line in Fig. 3 shows the effect of turning off the 27r-exchange 
three-body terms eqs. (6,7,13,14) in our calculation. After such a truncation there is almost no 
trace of nuclear matter saturation left over. In order to keep still a very shallow minimum at 
kfQ = 262 MeV the monopole mass has been reduced to A = 0.77 GeV. Leaving Aq = 1.15 GeV 
unchanged the curve for E{kf) without three-body terms would just decrease monotonically for 
all densities < p < 0.5 fm~^, reaching the huge negative value —85 MeV at its lower end. 

The pion-nucleon monopole form factor F{q^) introduced in our calculation is only a reg- 
ulator eliminating high momentum components. It is not a physical (observable) quantity. 
Interpreted differently, the monopole mass A is a control parameter which monitors the "spatial 
resolution" at the pion-baryon interaction vertices. A low-energy quantity, such as the nuclear 
matter equation of state, should not depend it. We take the point of view, characteristic of an 
effective field theory, that the explicit pion-exchange terms are accompanied by additional "un- 
resolved" short-distance contributions. The most general contact-interaction of four nucleons 
(momentum-independent and quadratic in momenta) gives rise to the following counterterm 
contribution to the energy per particle: 

k^ 

£(M'"" = ft(A)4 + B=(A)4. (17) 

with S3 (A) and -B5(A) two dimensionlcss coupling strengths. In the sense of a "renormalization 
group equation" we allow these contact-couplings to run with the resolution scale A such that the 
physical observable, the total energy per particle E{kf), is (as closely as possible) independent 
of A. 
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Figure 5: Energy per particle E{kf) of isospin-symmetric nuclear matter versus the nucleon 
density p = 2A;j/37r^. The curves correspond to three different monopole masses A = 0.5 GcV, 
1.0 GcV and 1.5 GeV with the contributions of the compensating contact-couplings ^3^5 (A) 
included. 



Fig. 4 shows how the two contact- couplings B^i^A) and -B5(A) must vary with the monopole 
mass A if the nuclear matter saturation point is required to stay fixed. The approximate linear 
rise of -B3(A) over the whole range 0.4 GeV < A < 2.4 GeV finds its explanation in the fact that 
for the iterated Ivr-exchange contributions eqs.(3,5), the monopole form factor serves also as a 
regulator on a linearly divergent loop integral (scaling as 1/r ~ A) which contributes to the 
energy per particle linear in density (p ~ /c^). 

In Fig. 5 we show the resulting nuclear matter equation of state for three different monopole 
masses A = 0.5 GeV, 1.0 GeV and 1.5 GeV. In the density region < p < 0.4 fm~^ there is almost 
no dependence on the resolution scale A left over. As expected, the small spreading of the three 
curves increases with increasing Fermi momentum kf, where higher momentum components of 
the interactions get probed more sensitively. It is also astonishing to observe that the rather 
complicated density dependence of the r-dependent terms (r = m.„/ A) written down in section 2 
can be almost perfectly counterbalanced by two leading powers of the Fermi momentum, and 
kf, reflecting contact-interactions without and with two derivatives. We have investigated the 
differences in the energy per particle with respect to the case Ao = 1.15 GeV where the running 
contact- couplings are zero: i?3(Ao) = -B5(Ao) = 0. For a wide range of monopole masses, 
0.4 GeV < A < 2.5 GeV, these deviations oscillate merely between -M.5MeV and -1.0 MeV in 
the whole density region < p < 0.5 fm~^. Therefore one can speak of a pretty stable nuclear 
matter saturation curve in our calculation which combines pion-exchange dynamics including a 
regularizing monopole form factor with (two) scale-dependent short-distance contact-couplings. 
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4 Chiral condensate at finite density 



The chiral condensate (0|gg|0) is an order parameter of spontaneous chiral symmetry breaking. 
With increasing temperature the chiral condensate decreases (or " melts away" ) . For low temper- 
atures this effect can be systematically calculated in chiral perturbation theory. At three-loop 
order [11] the estimate Tc ~ 190 MeV of the critical temperature, where chiral symmetry will be 
eventually restored, has been found. This (extrapolated) value of Tc is remarkably consistent 
with the recent result Tg = (192 ± 8) MeV obtained in numerical simulations of full QCD on the 
lattice [12]. The chiral condensate drops also with increasing baryon density. The leading linear 
term in the nucleon density p is readily derived (via the Feynman-Hellmann theorem) by differ- 
entiating the energy density of a nucleonic Fermi gas, pMjy + 0{p^^^), with respect to the light 
quark mass rriq. This introduces the nucleon sigma-tcrm ajy = {N\mqqq\N) = rnqdM^ / diriq as 
the driving term for the density evolution of the chiral condensate. Corrections beyond the lin- 
ear density approximation arise from the A'^A'^-interactions which transform the nucleonic Fermi 
gas into a nuclear Fermi liquid. Because of the Goldstone boson nature of the pions, ~ m^, 
the explicit pion-exchange dynamics in nuclear matter plays a particularly important role for 
the in-medium chiral condensate {qq){p)- Converting quark mass derivatives into pion mass 
derivative and using the Gell-Mann-Oakes-Renner relation m^fl = —rriq {0\qq\0) one finds for 
the ratio of the in-medium to vacuum chiral condensate: 



(99) (p) P 



{Q\qq\Q) 2mJl 

The nucleon sigma-term: 



2(Tn ^ dE{kf) 



(18) 



a ^M^^H^^Me. (19) 

* druq 2 

measures that portion of the nucleon mass Mjv which arises from the explicit chiral symmetry 
breaking in QCD (i.e. the non- vanishing up/down quark mass niq ^ 0). The empirical value of 
(77V as extracted from dispersion relation analyses of TrA^-scattering data is — (45 ± 8) MeV 
[13]. 

From our nuclear matter calculation with explicit one- and two-pion exchange dynamics we 
can now easily compute the derivative dE{kf)/dmT^ of the energy per particle with respect to 
the pion mass. In addition to the explicit m7r-dependence we consider also the implicit one 
through the nucleon mass: OMj^/dm-,, = 2aN/mn — 0.67. The smaller and less well-known 
quark mass derivatives of the axial-vector coupling constant qa and the pion-decay constant f^^ 
are neglected at present. For example, one-loop chiral perturbation theory gives for the pion- 
decay constant the small derivative: dfT^/dm.,, = {U — l)m^/(87r^/7r) — 0.063 [14]. Furthermore, 
from the combination of lattice QCD results for qa and their chiral extrapolation [15] one 
estimates 171^ dgA/drriTr < 0.05 around the physical point. Corrections to dE{kf)/dmTr from the 
mTT-dependence of the 7rA"-coupling constant ~ Qa/ fn, are therefore expected to be small, less 
than 5%. 

Fig. 6 shows the dropping chiral condensate (qq) (p) as a function of the nucleon density 
up to p = 0.4 fm"^. The full line corresponds to the linear density approximation using the 
central value ctat = 45 MeV. The other curves include the effects of nuclear correlations through 
the derivative dE{kf) / dniT^ for various values of the monopole mass A. We have assumed that 
the short-distance contact-couplings, B-i{A)/M'^ and B^{A)/M^, which keep the saturation 
curve E{kf) stable (see section 3) are quark mass independent. This is a natural assumption, 
but difficult to substantiate further. One observes in Fig. 6 that up to equilibrium density po 
the corrections beyond the linear density approximation are in fact very small. This finding 
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Figure 6: Density dependence of the chiral quark condensate {qq){p)- 

can be taken posteriori justification of the assumption made in ref.[16] about the in- 

medium behavior of the scalar mean-field S^'* . At higher densities a trend counteracting chiral 
restoration sets in. Similar features have been observed in the earlier chiral approach of Lutz et 
al. [1] as well as in ref.[17] where the relativistic scalar- vector mean-field phenomenology has been 
combined with estimates of the quark mass derivatives. It is comforting to see that the explicit 
one- and two-pion exchange dynamics does not lead to the opposite trend, namely acceleration 
of chiral restoration, which would undermine the foundation of the present approach to nuclear 
matter, based on the spontaneous breaking of chiral symmetry (in the vacuum). 

Finally, one should note that the nucleon sigma-term crjv has presently a sizable uncertainty 
of ±18%. Therefore, the error band associated with the linear density approximation 1 — 
P(7m/ (mlf^) masks practically all effects from nuclear correlations at least up to nuclear matter 
saturation density, po = 0.16 fm~^. 

5 Concluding remarks and outlook 

In this work we have performed a nuclear matter calculation which treats the long-range corre- 
lations from one- and two-pion exchange explicitly. In accordance with the recent approach of 
rcf.[5] based on the Viow-k potential we find that repulsive three-body terms are essential in order 
to achieve realistic binding and saturation of nuclear matter. As a novel feature we have intro- 
duced a pion-nucleon monopole form factor as a regulator to eliminate high-momentum compo- 
nents from the interactions. The dependence on the "resolution" parameter, the monopole mass 
A, can be perfectly counterbalanced by two running short-distance contact-couplings i?3 5(A). 
The resulting nuclear matter equation of state E{kf) is stable against variations of A. We have 
taken the pion mass derivative dE{kf) / duiT^ and obtained small corrections to the in-medium 
chiral condensate {qq){p) beyond the linear density approximation. 

A more detailed comparison of the present nuclear matter calculation with the one based on 
the Viow-k potential [5] offers prospects for a better understanding of the relevant short-range 
A^A^-dyn amies, which so far enters in the form of two adjusted contact-couplings S3 5(A). Work 
along these lines is in progress [18]. 
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